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ABSTRACT 

Recently, the Galactic Center has been reported to be a source of very high energy (VHE) 7-rays by the 
VERITAS, CANGAROO and HESS experiments. The energy spectra as measured by these experiments 
sho'w substantial differences. In this Letter 'we present MAGIC observations of the Galactic Center, result¬ 
ing in the detection of a differential 7-ray flux consistent 'with a steady, hard-slope po'wer law, described as 
dW.^/(dAdMi?) = (2.9 ± 0.6) x 10“^^(£’/TeV)“^'^®*^°'^ cm“^s“^TeV“^. The 7-ray source is centered at 
(Ra, Dec)=(17**45®®‘20®, -29°2'). This result confirms the previous measurements by the HESS experiment and 
indicates a steady source of TeV 7-rays. We briefly describe the observational technique used, the procedure 
implemented for the data analysis, and discuss the results in the perspective of different models proposed for 
the acceleration of the VHE 7-rays. 

Subject headings: gamma rays: observation — Galaxy: center — acceleration of particles 


1. INTRODUCTION 

The Galactic Center (GC) region contains many remark¬ 
able objects which may be responsible for high-energy 
processes generating Y -rays ilAharonian & NeronovI l2flfl-5t 
lAtovan & nermei1l2004ll . The GC is ric h in massive stel¬ 
lar cl usters with up to 100 OB stars JMorris & Serahviil 
[HU), immersed in a dense gas. There are young super¬ 
nova remnants, e.g. Sgr A East, and nonthermal radio arcs 
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JLaRosa et alJ200(il) . The dynamical center of the Milky Way 
is associated with the compact radi o source Sgr A*, which i s 
believed to be a massive black hole (iMorris & Serabvnl199H) . 
Within a radius o f 10 pc around the GC there is a mass of 
about 3 • IO^Mq JSchodel et alJ2002l . 

EGRET has detected a strong source in the direction o f 
the GC, 3 EG J1745-2852 l [Maver-]Hasselwander et alHT^ . 
which has a broken power law energy spectrum extending up 
to at least 10 GeV, with a spectral index of 1.3 below the break 
at a few GeV. Assuming a distance of the GC of 8.5 kpc, the 
7-ray luminosity of this source is very large, 2.2 • lO^^erg/s, 
which is equivalent to about 10 times the 7-flux from the Crab 
nebula. However, an inde pendent analysis of the EGRET data 
dHooper & Dingusll2002h indicates a point source whose po¬ 
sition is different from the GC at a co nfidence lev el beyond 
99.9 %. This was recently sustained bv gohil 12001 . 

In very high energy 7-rays the GC has been observed 
by VERITAS, CANGAROO and HESS (iKosack et aT]l2()()l 
iTsuchiva et al.l l2004t lAharonian et alJ 120041) . The energy 
spectra as measured by these experiments show substantial 
differences. This might be due to different sky integration re¬ 
gions of the signal or a source variability at a time-scale of 
about one year. 

2. OBSERVATIONS 

MAGIC (see e.g., iBaixeras et alJ d2004l) : ICortina et al.l 
l2()()l for a detailed description) is currently the largest single 
dish Imaging Air Cherenkov Telescope (lACT) in operation. 
Located on the Canary Island La Palma (28.8°N, 17.8°W, 
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TABLE 1 
Data Set 


Period 

date 

ZA n 

time [h] 

events [10®] 

obs. mode 

I 

Sep. 2004 

62-68 

2 

0.8 

ON 

II 

May 2005 

58-62 

7 

2.8 

wobble 

III 

lun./Jul. 2005 

58-62 

17/12 

6.4/5.0 

ON/OFF 


Note. — Data set per observation period of the GC. The column “time” states the 
effective observation time, the column “events” states the events after image cleaning 


2200 m a.s.l.), the telescope has a 17-m diameter tessellated 
parabolic mirror, supported by a light weight carbon fiber 
frame. It is equipped with a high quantum efficiency 576- 
pixel 3.5° field-of-view photomultiplier camera. The analog 
signals are transported via optical fibers to the trigger elec¬ 
tronics and are read out by a 300 MSamples/s FADC system. 

At La Palma, the GC ((RA, Dec) = (17‘'45“36", -28°56')) 
culminates at about 58° zenith angle (ZA). The star field 
around the GC is non-uniform. In the region west of the 
source (RA > RAgc + 4.7™) the star field is brighter. Within 
a distance of 1 ° from the GC there are no stars brighter than 
8 *^ magnitude. 

The MAGIC observations were carried out in the ON/OFF 
mode as well as in the false-source tracking (wobble) mode 
dFomin et al.lll994ll . The sky directions (Wl, W2) to be 
tracked in the wobble mode are chosen such that in the cam¬ 
era the star field relative to the source position (GC) is similar 
to the star field relative to the mirror source position (anti¬ 
source position); W1AV2 = (RAqo^ Decoc i 0.4°). During 
one wobble mode data taking, 50% of the data is taken at 
Wl and 50% at W2, switching between the two positions ev¬ 
ery 20 minutes. Dedicated OFF data have been taken, with 
a sky field similar to that of the ON region. The OFF re¬ 
gion is centered at the Galactic Plane, GCqff = (RA, Dec) = 
(17’°51™12", -26°52W'). In the same night OFF data was 
taken directly before or after the ON observations under the 
same weather conditions and hardware setup. 

After initial observations in September 2004 the GC was 
observed for a total of about 24 hours in the period May-July 
2005. Table^summarizes the data taken. 



Fig. 1.— Smoothed sky map of 7 -ray candidates (background subtracted) 
in the direction of the Galactic Center for SIZE > 300 ph. el. (corresponding 


to an energy threshold of about 1 TeV). Overlayed are green contours (0. 
Jy beam“^) of 90 cm VLA (BCD configuration) radio data ILaRosa et a 


) of 90 cm VLA (BCD configuration) 
The white line shows the galactic plane. 
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3. DATA ANALYSIS 

The data analysis has been carried out using the 
standard MAGIC ana lysis and reconstruction software 
dBretz & Wagnerll2003l) . th e first step of whi ch involves the 
calibration of the raw data (iGaug et alJl2005l) . After calibra¬ 
tion, image cleaning tail cuts of 10 ph. el. (core pixels ) and 5 
ph. el . (boundary pixels) have been applied (see e.g. (iFeganI 
These tail cuts are accordingly scaled for the larger 
size of the outer pixels of the MAGIC camera. The camer a 
images are parameterized by image parameters dFfillasl198.5h 

In th is analysis, the Random Forest method (see lBock et al.l 
d2004l) for a detailed description) was applied for the 7 /hadron 
separa t ion an d the energy estimation (for a review see e.g. 
iFe^ffW^ 'l. For the training of the Random Forest a sam¬ 
ple of Monte Carlo (MC) generated 7 -showers was used to¬ 
gether with about 1 % independent and identically distributed 
events drawn from the measured OFF-data. The MC 7 - 
showers were generated between 58° and 68 ° ZA with en¬ 
ergies between 10 GeV and 30 TeV. For the analysis of the 
September 2004 data set the Random Forest cuts were deter¬ 
mined using a sub-set of Galactic OFF data as background. 

The source-position independent image parameters SIZE, 


Fig. 2.— Distributions of 9 ^ values for the source and anti-source, see 
text, for SIZE > 300 ph. el. (corresponding to an energy threshold of about 
1 TeV). 

WIDTH, LENGTH, CONC llHillaslll985h and the third mo¬ 
ment of the ph. el. distribution along the major image axis, as 
well a s the source-position dependent parameter DIST llHillasI 
dUl), were selected to parameterize the shower images. Af¬ 
ter the training, the Random Eorest method allows to calcu¬ 
late for every event a parameter, called hadronness, which is 
a measure of the probability that the event is not 7 -like. The 
7 -sample is defined by selecting showers with a hadronness 
below a specified value. An independent sample of MC 7 - 
showers was used to determine the efficiency of the cuts. 

The analysis at high zenith angles was developed and veri¬ 
fied using Crab data with a ZA around 60°. The determined 
Crab energy spectrum was found to be consistent with other 
existing measurements (see Eig.|3 dot-dashed line). 

Eor each event the arrival direction of the prim ary in sky co¬ 
ordinates is estimated by using the DISP-method ilEomin et al J 
ll994llLessard et Tll200U iDomingo-Santamaria et al.ll200-5ft . 

Eor the sky map calculation only source independent image 
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parameters are used in the Random Forest training. Figure 
E] shows the sky map of o'-ray candidates (background sub¬ 
tracted, see e.g. llRowel]l2003l) ') from the GC region (obserya- 
tion periods II/III). It is folded with a two-dimensional Gaus¬ 
sian with a standard deyiation of 0 . 1 ° (roughly corresponding 
to the MAGIC PSF) and height one. A lower SIZE cut of 300 
ph. el. has been applied, corresponding to an energy threshold 
of about 1 TeV. The sky map is oyerlayed with contours (0.3 
Jy b eam~^l of 90 cm VL A (BCD configuration) radio data 
from iLaRosa et alJ(l 2000 l) . The brightest non-central source is 
the Arc. The excess is centered at (RA, Dec) = (17H5“20®, 
-29°2') (J2000 coordinates). The systematic pointing uncer¬ 
tainty is estimated to be 2 ' (for descr i ption of the MAGIC 
telescope driye system see iBretz et ^ ll2003h l and might in 
the f uture be further red uced with the MAGIC starfield mon¬ 
itor dRiegel et alJl200-'i^ . The excess is compatible with a 

S ioint source emission . The VHE 7 -ray source G 0.9H-0.1 
Aharonian et alj2005l) is located inside the MAGIC field-of- 
yiew. It sh ows a small excess con sistent with the low flux 
reported by lAharonian et akl d2005l) . The MAGIC excess is 
not yet statistically significant for the giyen exposure time. 

Eigure|2] shows the distribution of the squared angular dis¬ 
tance, 9^, between the reconstructed shower direction and the 
nominal GC position corresponding to Eig.n(obseryation pe¬ 
riods II/III). The obseryed excess in the direction of the GC 
has a significance of 7.3 standard deyiations (9^ < 0.02°). 
The source position an d the flux leyel are con sistent with the 
measurement of HESS ( lAharonian et alj2004ft within errors. 



Fig. 3.— Reconstructed VHE 7 -ray energy spectrum of the GC (statis¬ 
tical errors only). The full line shows the result of a power-law fit to the 
data points. The dash ed line shows the result of the HESS collaboration 
lAharonian et '^l 200 l . The dot-dashed line shows the ener gy spectrum of 
the Crab nebula as measured by MAGIC IWagner et all200'^ . 

Eor the determination of the energy spectrum, the Ran¬ 
dom Eorest was trained including the source dependent im¬ 
age parameter DIST with respect to the nominal excess posi¬ 
tion. Eor the spectrum determination only the largest data set 
(period III) was used. The cut on the hadronness parameter 
(50% 7 -efficiency corresponding to an effectiye area of about 
250000 m^) resulted in about 500 excess eyents with a min¬ 
imum SIZE of 200 ph. el. Eigure |3 shows the reconstructed 
VHE 7 -ray energy spectrum of the GC af tp the unfolding 
with th e instrumental energy resolution, see iMizobuchi et all 
( l2005l) . The differential 7 -flux can be well described by a 


simple power law: 

^j^ = (M±0.6)xlO-7£/TeV)--«- 

cm-7"^TeV"^ . 

The giyen errors (Icr) are purely statistical. The systematic 
error is estimated to be 35% in the flux leyel determination 
and 0.2 in the spectral index. 

Eigure |4] shows the reconstructed integral VHE 7 -ray flux 
aboye 1 TeV as a function of time. All GEE data are used 
for each time bin resulting in some correlation between the 
time bins. Different obseryation modi may result in differ¬ 
ent systematic errors. The flux leyel is steady within errors 
in the time-scales explored within these obseryations, as well 
as in the two year time-span between the MAGIC and HESS 
obseryations. 



time [Mjd] 


Fig. 4.— Light curve: Reconstructed integral VHE gamma-ray flux above 
1 TeV as a function of time. Within errors {Icr) the data are consistent with a 
steady emission. 

4. DISCUSSION 

Recent obseryations of TeV 7 -rays from the GC confirm 
that this is a yery important region for high energy processes 
in the Galaxy. In fact, this is not surprising since many dif¬ 
ferent objects, able to accelerate particles aboye TeV ener¬ 
gies, are expected there. The most likely source seems to be 
the massiye black hole identified with Sgr A* due to the di¬ 
rectional consistency. A blazar-like relatiyistic jet originating 
from the spinning GC black hole m ight be expected to pro¬ 
duce TeV 7 -rays (lEalcke et aklll 993l) . but flux predictions of 
this model are on the low side due to an unfayorable orienta¬ 
tion of the j et axis. Moreoyer, a short- term yariability would 
be expected. lAtoyan & Dermeil il2004ft propose that electrons 
can be accelerated to sufficiently high energies at the termi¬ 
nation shock of the sub-relatiyistic wind from the central part 
of the adyection dominated accretion flow onto the GC black 
hole, in analogy to the pulsar wind nebulae. The authors ex¬ 
plain the broad band emission from Sgr A* (from radio to TeV 
7 -rays) and suggest that the GeV source obseryed by EGRET 
has another origin. This is consistent with the recent determi- 
na tion of the position of the EGRETsource 3EG J1746-2851 
by lHooper & Dingul(l2002h and lPohll ( 120051) . Other scenarios 
for the 7 -ray production in the yicinity of Sgr A*, both lep- 
tonic and hadronic, haye also been found to be consistent with 
the TeV obseryations (for reas onable sets of parameters) bu t 
not with the GeV obseryations (lAharonian & Neronoyl200^ . 
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It is generally expected that 7 -rays produced in such compact 
source models should show relatively fast variability. The 
same level of TeV flux reported by HESS in 2004 and by 
MAGIC in 2005, and also during their own observation pe¬ 
riods extending over a few months, rather suggest a stable 
source on a year time scale. However, the 7 -ray flux above 
2.8 TeV (3.7cr significance) reported by Whipple during the 
extended period from 1995 through 2003 is a factor 2 larger 
dKosack et al 1200^ . 

The origin of 7 -ray emission in other types of sources is 
also possible as demonstrated by the detection of the sec¬ 
ond TeV 7 -ray source in the direction of the GC consistent 
with the lo cation of the cornposite s upernova rein nant SNR 
G 0.9H-0.1 llAharonian et al.l ^2005^ . iPohil lll997h proposed 
that the GeV emission can be related to the GC rad io arc. 
iCrocker et all ( 120051) . see also lFatuzzo & Melial d2003ft . argue 
for the GeV and TeV emission coming from different sites of 
the shell of the very powerful supernova remnant Sgr A East. 

More extended 7 -ray emission might also originate in the 
interaction of relativistic particles with the soft radiation and 
matter of the central stellar cluster around the GC. These 
particles can be accelerated by e.g. a very energetic pul¬ 
sar, a 7 -ray burst source, shocks in the wi nds of the mas¬ 
sive stars, or a shell type supernova remnant dBednarekl2002t 
Biermann et alJl20041 lOuataert & Loebll2005t ICrocker et alJ 
2005tlGrasso & Maccionell2005l) . If the TeV 7 -rays are pro¬ 
duced by leptons scattering off the infrared photons from 
the dust heated by the U V stellar radiation (as discussed by 
lOuataert & LoebI ll2005h '). then the 7 -ray power at ~ 100 
GeV should be almost an order of magnitude higher due 
to scattering of UV radiation. The 7 -ray energy spectrum 
should steepen between ~ 0.1 — 1 TeV. Instead, the HESS 
collaboration re ports a simple power l aw spectrum between 
~ 0.2 — 10 TeV ( lAharonian et al.l2004ft . In order to produce a 
gamma-ray spectrum well described by a single power law up 
to ^ 20 TeV, hadrons should have energies of about 10^ TeV. 
Such ha drons diffuse through t he region of the TeV source 
(< 7 pc ■ lAharonian et al l iEool i on a time scale of the order 
of 10 "*^ years, assuming the average magnetic held strength 
in this region of 10“^ G and the Bohm diffusion coefficient. 
Therefore, the natural source of relativistic hadrons seems to 
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be the supernova remnant Sgr A East or the energetic pul¬ 
sar created in the supernova explos ion JCrocker et all 1200 ^ 
iLaRosa et alJ200-^lBednareki2002h . However, this relatively 
young source of relativistic hadrons cannot be identified with 
the last 7 -ray bur st in the center of our Galaxy if it appeared 
^ 10® years ago JBiermann et al inoi. 

The GC can also be the brightest so urce of VHE 7 - 
rays from particle dark matter annihilation JPrada et al.l2004 
iHooner et alJl2004llFlixll2005l) . Most SUSY dark matter sce¬ 
narios lead to a cut-off in the 7 -ray energy spectrum below 
10 TeV. The observed 7 -ray energy spectrum extends up to 
20 TeV. Thus most probably the main part of the ob served 
7 -rad iation is not due to dark matter annihilation ( iHornsI 
l2004l) . However, an extended 7 -ray source due to dark mat- 
ter annihilation peaking in t he region 10 GeV to 100 GeV 
llElsasser & Mannheiml200^ cannot be ruled out yet. 

5. CONCLUDING REMARKS 

The MAGIC observations confirm the VHE 7 -ray source at 
the Galactic Center. Th e measured flux is comp atible with the 
measurement of HESS dAharonian et al.ll2flfl4t) within eiTors. 
The VHE 7 -ray emission does not show any significant time 
variability; our measurements rather affirm a steady emission 
of 7 -rays from the GC region. The excess is point like, it’s 
location is spatially consistent with SgrA* as well as SgrA 
East. 

The nature of the source of the VHE 7 -rays has not yet been 
identified. Future simultaneous observations with the present 
Cherenkov telescopes, the GLAST telescope and in the lower 
energies will provide much better information on the source 
localization and variability of emission. This will shed new 
light on the nature of the high energy processes in the GC. 
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